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1.  Introduction 

In  data  networks  with  bursty  users  whose  population  may  be  changing,  the  most  appropri¬ 
ate  multiple-access  techniques  are  those  belonging  to  the  class  of  Random  Access  Algorithms 
(RAAs).  The  latter  do  not  require  knowledge  of  the  user  population  in  their  operations,  they  can 
accommodate  bursty  independent  users  with  relatively  low  delays,  and  they  can  induce  satisfac¬ 
tory  throughput. 

The  important  performance  characteristics  of  some  RAA  are:  Stability,  throughput, 
transmission  delays,  and  insensitivity  to  channel  errors.  Stability  ensures  operation  in  the  pres¬ 
ence  of  changing  user  population,  good  throughput  translates  to  satisfactory  utilization  of  the 
channel  bandwidth,  and  low  delays  are  self  explanatory.  Insensitivity  to  channel  errors  means 
maintaining  good  performance  characteristics  when  operating  in  noisy  environments.  The 
source  of  noise  may  be  interferences  from  other  transmissions,  multi-path  fading,  etc.,  and  its 
effects  can  be  modelled  by  erroneously  read  feedbacks,  [4], 

In  this  paper,  we  consider  the  RAA  in  [1]  and  [3],  and  we  present  its  analysis  in  the  pres¬ 
ence  of  two  simultaneous  types  of  feedback  errors.  The  latter  algorithm  operates  in  both  full 
sensing  and  limited  sensing  feedback  environments,  it  is  synchronous,  it  requires  CNC  binary 
(collision  versus  noncollision)  feedback  per  slot,  and  in  the  absence  of  feedback  errors  and  the 
presence  of  the  limit  Poisson  user  model,  it  attains  throughput  0.43  (same  as  the  algorithm  in 
[2]).  This  algorithm  also  has  simple  operational  properties  which  allow  its  analysis  in  the  pres¬ 
ence  of  strict  delay  limitations  [1],  and  induces  delays  that  are  uniformly  better  than  those 
induced  by  the  algorithm  in  [2].  As  first  found  in  [3],  and  as  will  be  further  exhibited  in  this 
paper,  it  also  has  superior  error  resistance  qualities. 

The  paper  is  organized  as  follows:  In  Section  II,  we  present  the  system  model  and  the 
description  of  the  algorithm.  In  Section  III,  we  include  detailed  stability  analysis.  In  Section  IV, 
we  draw  conclusions  and  discuss  comparisons  with  other  existing  algorithms. 

n.  The  Model  and  the  Algorithm 

We  consider  independent,  identical,  and  packet  transmitting  users,  with  stationary  and 
memoryless  packet  generating  processes,  who  communicate  with  each  other  via  a  single  com¬ 
mon  channel.  We  also  consider  a  synchronous  system,  where  the  channel  time  is  divided  into 
slots,  each  of  length  identical  to  the  length  of  a  packet,  and  where  each  packet  transmission  can 
only  start  at  the  beginning  of  some  slot.  We  c^l  a  slot  empty  (E),  successful  (S),  or  collision 
(C),  if  it  is  not  occupied  by  some  packet  transmission(s),  or  is  occupied  by  a  single  packet 
transmission,  or  at  least  two  simultaneous  packet  transmissions  have  occured  in  it.  respectively. 
We  assume  that  in  the  event  of  a  collision,  all  the  involved  packets  are  destroyed,  and 
retransmission  is  then  necessary.  We  consider  the  case  where  feedback  per  slot  exists,  and  is 
received  by  the  users  at  the  end  of  each  slot  without  propagation  delays.  We  assume  that  the 
observed  by  the  users  feedback  is  CNC  binary  (collision  versus  noncollision),  where  an  observed 
NC  slot  may  be  either  empty  (E)  or  successful  (S).  We  also  assume  that  due  to  noisy  conditions, 
the  following  feedback  errors  exist:  An  E  slot  may  be  seen  by  the  users  as  a  C  slot,  with  proba¬ 
bility  e.  An  S  slot  may  be  seen  by  the  users  as  a  C  slot,  with  probability  6.  A  C  slot  is  always 
read  correctly  as  a  C  slot  by  all  the  users. 

For  the  system  described  above,  we  adopt  the  RAA  in  [3].  This  RAA  utilizes  a  window  A 
on  the  arrival  axis,  it  is  implemented  independently  by  each  user  in  the  system,  and  has  a  full 
sensing  and  a  limited  sensing  versions.  The  full  sensing  version  requires  that  each  user  know  the 


total  feedback  history  of  the  system,  at  all  times.  The  limited  sensing  version  only  requires  that 
each  user  observe  the  feedback  continuously,  from  the  time  he  generates  a  packet  to  the  time  this 
packet  is  successfully  transmitted.  For  completeness,  we  describe  here  both  versions  of  the  algo¬ 
rithm.  In  our  descriptions,  time  will  be  measured  in  slot  units,  where  slot  t  occupies  the  time 
interval  [t,  t-^l).  We  will  then  denote  by  Xt  the  feedback  of  slot  t,  as  seen  by  the  users  in  the  sys¬ 
tem,  where  Xt  is  either  C  or  NC. 


The  Full  Sensing  Algorithm 

The  algorithm  utilizes  a  window  of  length  A.  Let  t  be  a  time  instant  such  that,  for  some 
ti<t,  all  the  packet  arrivals  in  (0,  ti]  have  been  successfully  transmitted  by  the  algorithm  and 
there  is  no  information  regarding  the  arrival  interval  (tj,  t],  and  such  that  t  corresponds  to  the 
beginning  of  some  slot.  The  instant  t  is  then  called  Collision  Resolution  Point,  (CRP),  the 
arrival  interval  (0,  ti]  is  called  "resolved  interval",  and  the  interval  (tj,  t]  is  called  "the  lag  at  t". 

In  slot  t,  the  packet  arrivals  in  (ti,  t2  ^  min  (ti  +  A,t)]  attempt  transmission,  and  the  arrival  inter¬ 
val  (ti,  t2]  is  then  called  the  "examined  interval".  If  Xt  =  NC,  then  (ti,  t2]  contains  at  most  one 
packet  and  is  resolved  at  t.  If  Xt  =  C,  instead,  then  (ti ,  t2]  either  contains  at  least  two  packets  or 
it  contains  at  most  one  packet  but  erroneous  feedback  has  been  observed.  In  the  latter  case,  a 
collision  at  t  is  perceived  by  the  users,  and  its  resolution  starts  with  slot  t+1.  Until  the  perceived 
collision  at  t  is  resolved,  no  arrivals  in  (t2,  «’)  are  allowed  transmission.  The  time  period 
required  for  the  resolution  of  the  latter  perceived  collision  is  called  the  Collision  Resolution 
Interval,  (CRI).  During  some  CRI,  each  user  acts  independently  via  the  utilization  of  a  counter 
whose  value  at  time  t  is  denoted  rt.  The  counter  values  can  be  either  1  or  2,  and  they  are  updated 
and  utilized  according  to  the  rules  below. 

1.  The  user  transmits  in  slot  t,  if  and  only  if  rt  =  1.  A  packet  is  successfully  in  t,  if  and  only  if 
rt  =  1  and  Xt  =  NC. 

2.  The  counter  values  transition  in  time  as  follows: 

(a)  If  Xt_i  =  NC  and  rt_i  =  2,  then  r,  =  1 

(b)  If  Xt_i  =  C  and  q-i  =  2,  then  rt  =  2 

(c)  If  X[_i  =  C  and  ri_i  =  1,  then 

1,  with  probability  0.5 

2,  with  probability  0.5 


Remark  1.  From  the  operations  described  above,  it  is  not  hard  to  see  that  a  CRI  which  starts  with 
a  C  slot  ends  the  first  time  two  consecutive  NC  slots  occur.  Even  in  the  erroneous  feedback 
environment  considered  here,  upon  such  occurence,  the  u.sers  "know"  with  certainty  that  some 
CRI  has  ended. 

The  Limited  Sensing  Algorithm 

Let  some  user  generate  a  new  packet  within  the  time  interval  [t]  t]  I).  Then,  he  immedi¬ 
ately  starts  observing  the  feedback  sequence  {xi)i>,j,  starting  with  the  feedback  .  Let  us 
define  the  sequence  {ti},22.  as  follows:  t2  is  the  first  time  after  ti,  such  that  x,^  =x,j_i  =  NC. 
Then,  as  explained  above,  t2  corresponds  to  the  ending  slot  of  some  CRI,  and  from  t2  +  1  on,  the 
user  can  identify  the  ending  slots  of  CRIs  induced  by  the  algonthm  (where  some  of  the  latter 
CRIs  may  have  unity  length).  Each  tj  corresponds  to  the  ending  slot  of  a  CRI,  and  t,+i  is  the  first 


after  t;  observed  by  the  user  such  point.  At  t;,  the  user  updates  his  arrival  instant  to  the  value 
^  ti  +  (i~2)A;  we  call  the  sequence  {t^‘^)i>2,  updates.  Let  t^  be  such  that:  tke{ti}i>2, 

t^'^  <  ti  -  1  -  A;  V  i  <  k-1,  and  t^’‘^  >  t^  -  1  -  A.  Then,  in  slot  t^+l,  the  user  enters  a  CRI,  and 
transmits  his  packet  successfully  during  its  process.  The  operations  of  the  algorithm  during  the 
latter  CRI  are  exactly  as  described  under  its  full  sensing  version. 

Remark  2.  From  the  above  description,  it  is  clear  that  the  following  global  operations  are 
induced  by  the  limited  sensing  version  of  the  algorithm:  Let  T  be  a  slot  that  corresponds  to  the 
end  of  some  CRI.  Then,  in  slot  T+1,  all  the  users  with  current  updates  in  (t-A-1,  T— 1]  transmit. 
If  X7+1  =  NC,  then  the  CRI  which  started  with  slot  T+1  lasts  one  slot,  and  a  new  CRI  starts  with 
slot  T+2.  If  xj+i  =  C,  instead,  then  a  collision  is  perceived  by  the  users  whose  resolution  starts 
with  slot  T+2.  No  arrivals  that  did  not  participate  in  the  perceived  collision  at  T+1  are  transmit¬ 
ted,  until  the  latter  is  resolved.  During  the  collision  resolution  process,  the  users  operate  as  in 
items  1  and  2,  in  the  description  of  the  full  sensing  version  of  the  algorithm. 

III.  Stability  Analysis 

For  the  stability  analysis,  we  will  adopt  the  limit  Poisson  user  model  (infinitely  many  Ber¬ 
noulli  users).  Indeed,  as  proven  in  [6],  the  throughput  obtained  under  this  user  model,  is  a  lower 
bound  to  the  throughput  in  the  presence  of  any  number  of  independent  and  identical  users  whose 
packet  generating  processes  are  memoryless.  We  note  that  both  the  full  sensing  and  the  limited 
sensing  versions  of  the  algorithm  attain  identical  throughputs,  in  the  presence  of  the  limit  Pois¬ 
son  user  model;  the  delays  induced  by  the  latter,  however,  are  generally  longer  than  those 
induced  by  the  former.  Thus,  in  our  stability  analysis,  it  suffices  to  consider  the  full  sensing  ver¬ 
sion  of  the  algorithm  only. 

Consider  the  system  model  in  Section  II,  and  the  full  sensing  version  of  the  algorithm.  Let 
the  system  start  operating  at  time  zero,  and  let  then  {Cj)i>i  denote  the  sequence  of  lags  induced 
by  the  algorithm.  The  first  lag  corresponds  to  the  empty  slot  zero;  thus,  Cj  =  1.  In  addition, 
{Q  )i>i  is  a  Markov  chain  whose  state  space  is  at  most  countable.  To  see  that,  let  us  define: 

/u,d:  The  number  of  slots  needed  to  examine  an  interval  of  length  u,  given  that  the  lag 

equals  d. 

P(/ 1  u,d):  Given  that  the  interval  to  be  examined  has  length  u  and  the  lag  equals  d,  the  pro¬ 

bability  that  the  corresponding  collision  resolution  interval  has  length  /. 

From  the  operations  of  the  algorithm,  we  then  have: 

J/d,d  ,  if  Cj  =d<  A  ;  w.p  P(/ld,d) 

^jCi-A  +  /A,d  ,  ifCi=d>A;  w.p.P(/IA,d) 

We,  therefore,  conclude  that  {Ci}i>i  is  a  Markov  chain  with  state  space  the  countable  set 
E  =  {x>l  :  X  =  k-mA,  k,  m  e  N}.  It  can  be  seen  that  any  state  can  be  reached  from  any  other; 
therefore,  {Ci}i>i  is  an  irreducible  Markov  chain.  Since  P(Ci+i  =  1  I  Cj  =  1)>0,  we  conclude 
that  {Cj}i>i  is  also  aperiodic.  Thus,  Pake’s  Lemma  [5]  applies,  and  gives  that  the  following 
condition  is  sufficient  for  the  ergodicity  of  the  Markov  chain: 


On  the  other  hand,  if  E(/IA,d)>  A,  then  the  Markov  chain  is  not  ergodic,  and  the  system  is 
unstable.  Let  denote  the  expected  length  of  a  collision  resolution  interval  given  that  it  starts 
with  a  collision  of  multiplicity  k.  We  can  then  write: 

E(/lA,d)=  X  LkC-^ (3) 

k=0 


In  Section  III.A,  we  show  that: 

(i)  Lie  ;  1^0  can  be  computed  recursively,  and 

(ii)  Lie  quadratically  upper  bounded.  Lie  <  LjJ  ^  a  k^  +  P  k+y ;  k>2. 

Due  to  the  upper  bound  on  L^,  we  conclude  then  that  the  following  condition  is  sufficient  for 
stability: 

30 

X  Lie  p(klA)+  X  Lk  p(klA)  <  A  14) 

k=0  k=31 


where 


p(klA)  =  e-^.S^ 

k! 


After  some  manipulations,  we  conclude  that  (4)  is  equivalent  to: 

30  f  30  ^  1 

f(XA)=  X  Lie  p(klA)  +  a-^  (XA)2  +  XA-  X  k  p(klA)^  + 


30  30 

+  p  XA- X  k  p(klA)  +Y  l-Xp(klA)  <A 
L  k=0  k=0 


Let  us  now  define: 


X  ^  XA 


Then,  from  (5)  and  (6),  we  conclude  that,  for  the  stability  of  the  algorithm,  it  is  sufficient  that  the 
input  rate  X  satisfies  the  following  inequality. 

X  <  sup  — ^  (7) 

f(x) 

The  following  condition  specifies  a  region  of  X  values  for  which  the  algorithm  is  unstable. 

X>sup  (8) 

g(x) 

30 

where  g(x)  =  X  e“’‘  x'‘/k! 
k=o 

The  maximization  of  expessions  in  (7)  and  (8)  has  been  done  numerically,  and  provides  the 
throughput,  as  well  as  the  optimal  window  size  A. 


III. A  Computation  of  Lv 
We  define: 

Gn,k-n-"  The  expected  number  of  slots  needed  by  the  algorithm,  for  the  successful  transmission 
of  k  packets,  given  that  n  of  those  packets  have  counter  values  equal  to  one,  and  k-n  of  the  pack¬ 
ets  have  counter  values  equal  to  two. 

Notice  that  =Gk,o.  for  k>2,  while  for  k<2.  We  first  show  how  we  compute  Lq 

and  L] . 

(i)  Computation  of  Lo: 

From  the  operation  of  the  algorithm  we  have 

1  ;  w.p.  (1-e) 


where 


From  (9a)  and  (9b)  we  find  that 


^  j  1  -f-  Go.o  ;  w.p.  e 

T  +Lo  ;  w.p.  (1-E) 
Go,o  -  j  I  +Go.o  ;  w.p.  e 


Lo  = 


1 


(l-e)"^ 


(ii)  Computation  of  Li 

It  was  found  that  Lj  satisfies  the  following 

1  ;  w.p.  (1-6) 

I 

Li  = 


where  Gi_o  andGo.i  satisfy  the  following 


Gi.o  - 


Go.i  =' 


From  (10a),  (10b),  and  (10c)  we  find  that 


1+Gi.o  i  w.p.  0.55 
l-(-Go.i  i  w.p.  0.55 

1-t-Lo  ;  w.p.  (1-6) 
I+Gi_o  ;  w.p.  0.55 
1-t-Go.i  i  w.p.  0.55 

1-t-Li  ;  w.p.  (1-E) 
I+Go,i  i  w.p.  E 


T  ,  — 

1  5  ■ 

-1 
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’  2-6 

l-t- 


2(1-e)  (2-5) 


(1-e)2  2(1-e)J 
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(10b) 
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(iii)  Computation  of  Lv,  for  k>2. 

From  the  operation  of  the  algorithm  we  obtain: 


where 


G„,k-„  =  1  +  S  i  2'"  Gi.k-i  ;  n>2  .  k>n 

i=0  ^ 


Gi,k  = 


l+Lk(l-S)  +  0.55  -^  +  Lk^i 
1-0.55 


It  can  be  shown  by  induction  that  Gn^k-n  has  the  following  form 

G„.k-n  =  A(i>Gk,o+Af  Gk-i.o  +  A(?)  ;  2^<k 

;  where  A^'\  i=l,2,3  are  independent  of  k  and  can  be  computed  recursively  as  follows: 

A  0)  -  2+5  *  n)  _  4(1-5)  H)  _  14-35-1 2e  +  4£5 

6-35  ’  6-35  ’  ^  3(l-e)  (2-5) 

Ay)  =  [l-2-"r'2-"|l+-^  +  "£[j]Apj-  ,  n>3 


Al('^  =  [l-2'"r‘  2' 


-m-i  I  n(l-5)  . 


(1-0.55)  ^2 


2:[?]APLn>3 


.  (3)  ^  r  1  o-n  1-1  J  1  .  — ^  ^  4. _ 5n^ 


For  n=k,  expression  (12)  gives: 

aP  a9^ 

'  i::XF 

Expression  (14)  together  with  the  recursions  in  (13)  provide  a  mean  for  the  computation  of 
Lk,  k>2. 

Development  of  an  upper  bou nd  on  Li. 


(lib) 


(13b) 


It  can  be  seen  by  induction  that: 


l-Af) 


=  1  ;  k>2 


(15b) 
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aP)  <  (2-2£-i-8)  1  .  . 

^  -  (2-5)(l-e)  (1-e) 

Here,  we  prove  (15b)  ;  the  proofs  for  (15a)  and  (15c)  are  similar.  By  summing  (13b)  and  (13c) 
and  after  some  simple  manipulations  we  obtain: 


+  A?)  = 


l+n+  X  I  -  I  (Ai*^  +  AP) 


2"  -1 


We  also  have  that  +  aP  =  1.  Next  we  assume  that  A['^  +  Ap^  =  1  for  i  =  3,4,..,  n-1. 
Expression  (16)  then  gives 


1+n  + 


A<»  +  A<?'  = 


n-l  r  ' 

2  i 

1—'^  ^ 


2"-l 


which  completes  the  proof. 


From  (14)  and  (15)  we  finally  find: 


,  3(2-2e+S)  3(2-5) 

^  (l-e)(4-55)  ^  (l-e)(4-55)  ’ 


from  which  we  can  show  that 


1  ^  3(2-2£+8)  2  ,  3(6-2e-5)  f 

2(l-e)(4-58)  2(l-e)(4-56)  * 


— ^  LlJ  =  otk^+pk+y  (19) 
(l-e)(4-55)  =  “  H  t 


III.B  Computation  of  Lv,  for  the  Capetanakis  Dynamic  Algorithm 

For  the  operation  of  the  algorithm  the  interested  reader  is  referred  to  [2]  and  [5].  Here,  the  quan- 
titites  Lk,  k>2  can  be  computed  recursively  as  follows: 


Lk  =  [l-2  2-''r‘-^  l+2Lo  2-'^  +  2 


k-l  fv-) 

XL.  t 

i=i  ^  '' 


where 


I  1  l-2e  +  5 

^  (l-2e)  (l-2e)(l-5)  ^ 

Moreover  the  following  upper  bound  on  Lk,k>4,  has  been  found. 

L  <  k  -  1  ^  Lj^  =  ^tk+v  (22) 

I  (l-2e)  (l-2e)(l-6)J  =  ^  ^  ’ 

From  (2),  and  due  to  the  upper  bound  in  (22),  we  conclude  that  the  following  condition  is  suffi¬ 
cient  for  stability: 

30  [30  1  r  30 

h(XA)=  X  Lk  P(klA)  +  ll  XA-Xkp(klA)  +v  l-Xkp(klA)  <A  (23) 

k=0  k=0  k=0 


or  equivalently 


.where  x  =  XA. 


[i] 

\u 


The  following  condition  specifies  a  region  of  the  input  rate  values  X,  for  which  the  algorithm  is 
unstable. 


30 

where  p(x)  =  X  P(k  I  A) 
k=o 

rV,  Conclusions  and  Comparisons 

We  analyzed  the  proposed  algorithm,  as  well  as  the  Capetanakis’s  Dynamic  Algorithm  [2], 
(in  [5],  the  Capetanakis’s  Nondynamic  algorithm  is  considered).  Binary  CNC  observed  feedback 
was  assumed.  Both  algorithms  attain  the  same  throughput  (X*  =  0.429)  when  e  =  5  =  0. 

In  Table  1,  the  throughputs  for  the  proposed  algorithm  and  the  Capetanakis’s  Dynamic 
Algorithm  are  presented  for  various  values  of  the  error  probabilities  e  and  8.  The  values  of 
e  and  5  have  been  chosen  to  demonstrate  that  the  proposed  algorithm  is  remarkably  insensitive 
to  feedback  channel  errors.  Even  for  the  practically  extreme  values  e  =  5  =  0.1,  the  throughput  is 
almost  90%  of  its  value  in  the  error  free  case.  From  the  results  presented  in  Table  1,  we  also 
conclude  that  the  proposed  algorithm  allows  operation  (positive  throughput)  as  long  as 
e  <  1  and  8  <  1,  while  if  e  >  0.5  the  throughput  for  the  Capetanakis’s  algorithm  is  then  zero.  We 
notice  for  example  that  for  e  =  0.5  and  8  =  0,  the  proposed  algorithm  attains  throughput  as  high 
as  X*  =  0.325. 
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